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We consider the dimensional crossover in the topological matter, which involves the transformation of dif- 
ferent types of topologically protected zeroes in the fermionic spectrum. In the considered case, the multiple 
Dirac (Fermi) point in quasi 2-dimensional system evolves into the flat band on the surface of the 3-dimensional 
system when the number of atomic layers increases. This is accompanied by formation of the spiral nodal lines 
in the bulk. We also discuss the topological quantum phase transition at which the surface flat band shrinks 
and changes its chirality, while the nodal spiral changes its helicity. 



1. INTRODUCTION 



Topological matter is characterized by nontrivial 
topology of the Green's function in momentum space 
[TJ [2]. The topological objects in momentum space 
(zeroes in the spectrum of fermionic quasiparticles) in 
many respects are similar to the topological defects in 
real space, and are also described by different homotopy 
groups including the relative homotopy groups. In par- 
ticular, the Fermi surface is the momentum-space ana- 
log of the vortex loop in superfluids/superconductors; 
the Fermi point (or Dirac point) corresponds to the 
real-space point defects, such as hedgehog (monopole) 
in ferromagnets; the fully gapped topological matter is 
characterized by skyrmions in momentum space, which 
are analogs of non-singular objects - textures; etc. 

Here we discuss the transformations of the topolog- 
ically protected zeroes, which occur during the dimen- 
sional crossover from a 2-dimensional to a 3-dimensional 
system. We consider the dimensional crossover which 
involves such topological objects as a nodal line in a 3D 
system; the flat band, which is an analog of the vor- 
tex sheet; and the topologically protected Dirac points 
■with multiple topological charge |iV| > 1 in quasi 2- 
dimensional substance, which are analogous to the mul- 
tiply quantized vortex. 

The Fermi bands, where the energy vanishes in a fi- 
nite region of the momentum space, and thus zeroes in 
the fermionic spectrum have co-dimension 0, have been 
discussed in different systems. The flat band appears 
in the so-called fermionic condensate [3l 01 [5j [6] . Topo- 
logically protected flat band exists in the spectrum of 



fermion zero modes localized in the core of some vortices 
[71|5J|5]- In particle physics, the Fermi band (called the 
Fermi ball) appears in a 2+1 dimensional nonrelativis- 
tic quantum field theory which is dual to a gravitational 
theory in the anti-de Sitter background with a charged 
black hole 10 . The flat band has also been discussed on 
the surface of the multi-layered graphene [TT] and on the 
surface of superconductors without inversion symmetry 
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The topologically protected 2-dimensional and 3- 
dimensional Dirac points with multiple topological 
charge N were considered both in condensed matter 
[13 HI E M US Hi El US El HQ] and for relativis- 
ts quantum vacua [5TJ [2 r23[ • In the vicinity of the 
multiple Dirac point with topological charge N the spec- 
trum may have the form E 2 cx p 2N . We consider the 
special model of the multilayered system discussed in 
|20) . where the topological charge N of the Dirac point 
coincides with the number of layers. In this model, when 
N — > oo, the multiple Dirac point transforms to the flat 
band in the finite region of the two-dimensional momen- 
tum on the surface of the sample. The interior layers 
in the limit N — > oo transform to the bulk state, which 
represents a semi-metal in which the nodal line (line 
of zeroes) forms a spiral. The projection of this spiral 
onto the edge layer produces the boundary of the flat 
band. The latter is similar to what occurs in supercon- 
ductors without inversion symmetry, where the region 
of the flat band on the surface is also determined by the 
projection of the topological nodal line in the bulk on 
the corresponding surface |12j . 
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2. FLAT BAND AND SPIRAL NODAL LINE 

Let us first consider the model in Ref. [5D] , specified 
in Sec. Q> .below, in the limit iV — > oo. The effective 
Hamiltonian in the 3-dimensional bulk system which 
emerges in the continuous limit N — > oo is the 2x2 
matrix 

H= (t o) ' f = P*- i Pv- t + e ~ iaP *- t - eiaPz - (!) 

Here t + = |< + |e i<#, + and t_ = |i_|e^~ are the hopping 
matrix elements between the layers and a is the inter- 
layer distance. The hopping matrix element propor- 
tional to a + (' is )■ The energy spectrum of the 
bulk system 

E 2 = [p x - \t+\ cos(ap z - 0+) - \t-\ cos(ap z - 0_)] 2 
+ [p y + \t+\ sm(ap z - <ft + ) - \t-\ sin(ap z - 0_)] 2 , 

(2) 

has zeroes on the line (see Fig. [TJ: 

p x = \t+\ cos(ap z - <j)+) + \t-\ cos(ap z - tj>J) , 
p y = \t-\ sm(ap z - <f>J) - \t + \ sin(ap z - (j) + ) . 

These zeroes are topologically protected by the topolog- 
ical invariant [5] 

2V"i = --^tr {cUa^H-^iH, (4) 

where the integral is along the loop around the nodal 
line in momentum space. The winding number around 
the element of the nodal line is N% = 1. For the inter- 
acting system the Hamiltonian matrix must be substi- 
tuted by the inverse Green's function at zero frequency, 
H — > G _1 (w = 0,p), which plays the role of effective 
Hamiltonian, see also |23j . 

The same invariant can be written if the contour of 
integration is chosen parallel to p z , i.e. at fixed pj_. Due 
to periodic boundary conditions, the points p z = ±n/a 
are equivalent and the contour of integrations forms the 
closed loop. 

JVi(pj_) = tr / + ' d Pz azH-^^H, (5) 

For interacting systems, this invariant can be repre- 
sented in terms of the Green's function expressed via 
the 3D vector g(p z ,w) [5]: 

G^{uJ,p z ) = ig z (uj,p z ) - g x (u,p z )cr x + g v (u>,p z )a y . 

(6) 




Fig. 1. Fermi line for the case t = |t+| = 10|i_| (red 
with circles) and t = = 10|t+| (blue) along with 
their projection to the p z = plane (dashed). This 
projection represents the boundary of the dispersionless 
flat band on the surface. In both cases (/>+ =0- — 0. 
Note that the helicity of the two lines is opposite, which 
gives the opposite signs of the invariant Ni(p±) = ±1 
in (J5]), and the opposite chiralities of the flat band. 

In our model the components g(p z ,w) are: 

g x (p z ,u>) = p x - \t+\ cos(ap z - 0+) - |t_| cos(ap z - 0_) 

g v (Pz,u) =p v + \t+\ sm(ap z - (j>+) - \t-\ sm(ap z - </>_) , 

g z (pz,u) = u . 

(7) 

Then the invariant §5§ becomes [5] 




where g = g/|g|. It describes the topological properties 
of the fully gapped ID system, with p x and p v being the 
parameters of the system. 

Let us first consider the case i_ = 0. For i_ = 0, 
the nodal line in ([3]) forms a spiral, the projection of 
this spiral on the plane p z = const being the circle 

= p x + p 2 = \t + \ 2 (the spiral survives for i_ ^ 0, 
but circle transforms to an ellipse, see Fig. [T] for the 
case |t_| < The topological charge in (J5]) is 

■^i(Px) = 1 for momenta |pj_| < If the momen- 

tum p^ is considered as a parameter of the ID system, 
then for |pj_| < \t+\ the system represents the ID topo- 
logical insulators. For |pj_| > \t + \ one has Ni(p±) = 
and thus the non-topological ID insulator. The line 
|pj_| = \t+\ marks the topological quantum phase tran- 
sition between the topological and non-topological ID 
insulators. 

Topological invariant iVi(px) in ([5]) determines also 
the property of the surface bound states of the ID sys- 
tem: the topological insulator must have the surface 
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Fig. 2. The Fermi line (nodal line) in the bulk at 
the topological quantum phase transition, which occurs 
at t — \t+\ = \t-\, for two values of the phase shift 
cf>- — 4>+. In this case the projection of the nodal line 
to the p z = plane shrinks to the line segment with 
zero area because the Fermi line is flat: it is within the 
corresponding plane drawn on the figure. As a result 
the dispersionless flat band on the surface is absent at 
the transition, and has opposite chiralities on two sides 
of the transition. 

states with exactly zero energy. These states exist for 
any parameter within the circle |pj_| = \t+\. This means 
that there is a flat band of states with exactly zero en- 
ergy, i?(|pj_| < = 0, which is protected by topol- 
ogy. The bound states on the surface of the system can 
be obtained directly from the Hamiltonian: 

H = <T x (p x - \t+\ cos(ap 2 )) + a y (p y + \t+\ sm(ap z )) 
p z = ~id z , z < . 

(9) 

We assumed that the system occupies the half-space 
z < with the boundary at z — 0, and made rotation 
in (p x ,Py) plane to remove the phase 0+ of the hopping 
element t+ . This Hamiltonian has the bound state with 
exactly zero energy, E(p±) = 0, for any |pj_| < \t + \, 
with the eigenfunction concentrated near the surface: 

* oc ) ( Px -i Py )e, P Z ^ P ; +ip ^ , | P J < \ U 

(10) 

The normalizable wave functions with zero energy exist 
only for p± within the circle |pj_| < i.e. the sur- 
face flat band is bounded by the projection of the nodal 
spiral onto the surface. Such correspondence between 
the flat band on the surface and lines of zeroes in the 
bulk has been also found in Ref. [T2"] . 

Restoring the non-zero hopping element we find 
that for \t-\ < \t+\ there is still the region of the mo- 



mentum p± for which the topological charge in (|5|) is 
Nx(p±) = 1. However, the area of the projection of 
the nodal line on the surface (and thus the area of the 
flat band) is reduced. Finally at |t_| — \t+\, the nodal 
line becomes flat, its projection on the surface shrinks 
to the line segment p y /p x = tan(0 + — 0_)/2, and the 
flat band disappears (Fig. [5]). For |t_| > \t + \, the spiral 
appears again, but the helicity of the spiral changes sign 
together with the topological charge in ([5]), which be- 
comes Ni(p±) = — 1. Thus the point = \t + \ marks 
the topological quantum phase transition, at which the 
flat band changes its orientation (or actually its chiral- 
ity). At the transition line \t-\ = the flat band 
on the surface does not exist. The non-zero helicity of 
the nodal line at |f_| ^ |f+| reflects the broken inver- 
sion symmetry of the system at |t_| 7^ |t+|. Note that in 
Ref. jH] the flat surface bands also appeared in systems 
without inversion symmetry. 



3. FORMATION OF THE FLAT BAND IN 
MULTILAYERED SYSTEM 

We consider the discrete model with finite number 
N of layers. It is described by the 2N x 2N Hamilto- 
nian with the nearest neighbor interaction between the 
layers in the form: 



Hij(p±) = 

a ■ pxSij - (t+<? + + t-a~)Sij + i - (t* + a~ 
1 < i < N , p ± = (p x , Py ) . 



+ t*_a + )S itj . 



(11) 



In the continuous limit of infinite number of layers 
transforms to (TT|). For t_ = and t+ = t, equation ([TT]) 
represents the particular case of the model [5^, which 
exhibits the Dirac point with multiple topological charge 
equal to N. 

Let us consider how this multiple Dirac point trans- 
forms to the flat band in the dimensional crossover, i.e. 
in limit N — > 00 (see Fig. [3]). The Hamiltonian H in 
(TIT]) with t- = and t+ = t has two low-energy eigen- 
'states with dispersion in the vicinity of multiple Dirac 
point (at |p x | < \t\) 



e/|t|«±(|p x |/|t|) 



(12) 



Let us look for the eigenfunctions corresponding to the 
dispersion (fT2")) . At p± = 0, the eigenfunctions are finite 
only in the first and TV'th layer, 



ip ± = ip(p± = 0) 



V2 



[U)i±|t)iv] 



(13) 
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Fig. 3. Formation of the surface fiat band for param- 
eters t = t+ G R, t- — 0. When the number N of 
layers increases, the dispersionless band evolves from 
the gapless branch of the spectrum, which has the form 
E — ±|px| in the vicinity of multiple Dirac point. 
The spectrum is shown as a function of p x for p y = 
and finite value of p z . The curves for N = 100 and 
N = 200 are almost on top of each other. Asymp- 
totically the spectrum E = ilpxl^ transforms to the 
dispersionless band within the projection of the nodal 
line to the surface. 



That is: we get Hipo± — 0. Now for |pj_| <C \t\ it suffices 
to find a function Sip satisfying HSip — \p±\ N /{t^^ 1 . 
It turns out that the small parameter in this expansion 
is dependent on the layer index, and the smallest cor- 
rection we have to include is rj ~ |p_i_| . In this order, 
we can expand the eigenvalue equation in rj: 



T] 2 Slp2 



Hip p = H(ip + vHi 
= rj{HSipi + r)HSip 2 + ■ 
= arjipo + ri 2 HSip2 = |p 
= V (Tjj /\t\ N - 1 +r,SiPi/\t\ N 



N /\t\ N -^ P 



(14) 



where r/Sipi = Sip. To lowest order in 77 it is thus suf- 
ficient to find this Sip. Moreover, to preserve the nor- 
malization, we seek Sip such that it is orthogonal to ipQ. 
In this case, to the lowest order in 77 the normalization 
remains unaltered. 

Let us hence write Sip in the form 



Sip 



N 

71=1 



I); 



(15) 



so that aii = a Nf — 0- Now acting with the Hamilto- 
nian yields 



N 



H5iP=Y. 



n=l 



(Px ~ iPy)a n t\ i)n + (Px + iPy)Gtnl\ t)r 



= ± 



1.5 



(pI+pI) n/ 



V2\t\ 



N- 



— (U>i ± I t>iv)- 



(16) 
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Fig. 4. Formation of the nodal line from the evolution 
of the gapped branch of the spectrum of the multilay- 
ered system, when the number TV of layers increases. 
The spectrum is shown as function of p x for p y = and 
finite value of p z for t = t+ = 2i_ £ R. The curves for 
N — 100 and N = 200 lie almost on top of each other, 
indicating the bulk limit. Asymptotically the nodal line 
Px = ±(t+ + t-) cos(ap z ), p y = ±(t+ -i_)sin(ap z ) is 
formed (two points on this line are shown, which corre- 
spond to p y = 0). 

Considering this equation separately for each com- 
ponent results in 



Or4 



± 



(-l)^- 1 ^ + l P y) N/2 {Px ~ l P y) N/2 - n 

V2/j( Ar - 1 )/ 2 (t*)( Ar + 1 )/ 2 -™ 

(-l)^- 1 ^ + iPy) n - N/2 (Px - *Py) N/2 
y /2 t n-(N-l)/2(t*)(N-l)/2 

The expressions are somewhat simpler for p y = 



Oinf = ± 



(-1) 



JV-l 



Px\ N - n ft 



(-1) 



V2 

JV-l 



(N-l)/2 



V2 



Px\ n ( t 



(N-l)/2 



(17a) 
(17b) 

(18a) 
(18b) 
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Wc hence find that the eigensolutions behave as ~ |p_i_| n 
for n layers away from the surfaces, in agreement with 
the wave function (fTU)) for the surface flat band obtained 
in the continuous limit 




Fig. 5. Gapped spectrum as a function of p x and 
p y for N=200 layers corresponding to the bulk limit. 
The nodal line forms at the projection of the Fermi line 
Px = ±(t+ + t-)cos(ap z ), p y = ±(t+ - t-)sin(ap z ) to 
the p z -plane. 



We have considered the evolution of the multiple 
Dirac point with spectrum E ~ ilpjj-^ into the surface 
flat band when N — > oo. However, the surface flat band 
survives when the coupling t- is added which splits the 
multiple Dirac point. The reason for the robustness of 
the flat band is the nodal line which is developed in 
the bulk (see Figs. 0] and [5]). As we discuss in the pre- 
vious section, the nodal line supports the topological 
stability of the flat band due to the bulk-surface cor- 
respondence. Formation of the dispersionless flat band 
for finite |t_| < \t + \ is shown in Figure HO The flat band 
has the ellipse-shaped region (Fig. [7]) whose boundary 
is the projection of the bulk nodal line on the surface. 




Px/t 



Fig. 6. Evolution of the spectrum at t = t+ = 2t- £ R. 
The curves for N = 100 and N = 200 are almost on 
top of each other, indicating the bulk limit. Asymptot- 
ically the flat band is formed for \p x \ < \t+\ + \t-\ and 
\Pv\ < \U\ - \t-\. 




Fig. 7. Energy of the surface states for different p x 
and p y for N = 200 layers corresponding close to the 
bulk limit, calculated with the same parameters as in 
Fig. [6] The flat band forms in the ellipse-shaped re- 
gion \p x \ < (t+ + t_)cos(0), \p y \ < (t+ - t-)sin(0), 
9 G [0,2tt]. 

4. CONCLUSION 

We have considered the dimensional crossover in 
which the multiple Dirac (Fermi) point in quasi 2+1 
system evolves into a flat band on the surface of the 
3+1 system when the number of atomic layers increases. 
The formation of the surface flat band is a generic phe- 
nomenon which accompanies the formation of the nodal 
lines in bulk in the form of a spiral. We have also demon- 
strated a new type of a topological quantum phase tran- 
sition, at which the flat band shrinks and changes its chi- 
rality. This transition is accompanied by the change of 
the helicity of the nodal line. The considered crossover 
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is one of the numerous examples of the evolution of the 
topologically non-trivial quantum vacua, which are rep- 
resented by the momentum-space topological objects. 
This example displays the ambivalent role of symmetry 
for these objects: the symmetry may support the topo- 
logical charge of the object in momentum space, or may 
destroy the object. In our case the time reversal sym- 
metry supports the existence of the nodal line in bulk, 
and thus the flat band on the surface, while the space 
inversion symmetry kills the flat band. 

Different scenarios of the dimensional crossover can 
be realized with cold atoms in optical lattices [24]. For 
our scenario we need the special stacking of graphene- 
like layers to have a spiraling nodal line. The spiral 
formed by zeroes in the energy spectrum has been dis- 
cussed in Ref. [25] for rhombohedral graphite. The 
nodal line there is modified to the chain of the con- 
nected electron and hole Fermi pockets, so that "the 
Fermi surfaces resemble very long 'link sausages' wound 
in a loose spiral" [2"5] . 

The dispersionless flat band also exists on the surface 
of the polar phase of triplet superfluid/superconductor. 
This superconductor obeys the time reversal and space 
inversion symmetry, and it has a line of zeroes in the 
form of a ring pQ. This ring gives rise to two surface 
flat bands with opposite chirality corresponding to two 
directions of spin. However, spin-orbit interaction may 
lead to the mutual annihilation of the flat bands. 

More examples of dimensional crossover and exotic 
quantum phase transitions emerge when one considers 
the topology in the phase space, i.e. in the combined 
momentum- real (p,r) space [261 IT]. This is appropriate 
in particular for the fermion zero modes localized on 
topological defects, which also may have dispersionless 
flat band [JJ UJ and bulk- vortex correspondence [37J [35] • 
The flat band inside the vortex core in 3-dimentional 
supcrfluids with Fermi (Dirac) points emerging due to 
the bulk- vortex correspondence is discussed in |29) . 

Systems with topologically protected Fermi lines or 
Fermi points belong to the broad class of topological 
matter. As distinct from topological insulators and su- 
perconductors/superfluids of the 3 He-B type [3U1 I5"T] . 
which belong to fully gapped topological matter, these 
are the gapless topological matter. However, it has the 
features which was earlier ascribed only to topological 
insulators, i.e. protected gapless states on the surface 
or inside the vortex core. If one or two components 
of the momentum p is fixed, such as the projection p z 
of the momentum on the direction of the vortex axis 
(see accompanied paper [25]) or \p±\ < \t\ in our case, 
the system effectively behaves as one-dimensional and 
two-dimensional topological insulator correspondingly. 



This is because for these parameters the system is fully 
gapped, while the effective ID or 2D Hamiltonian has 
a non-trivial topology. Since the topological insulators 
cannot be adiabatically turned to a trivial insulator, this 
gives rise to the zero energy edge states in those inter- 
vals of parameters p z or \p± | , for which the topology is 
nontrivial. As a result, in both cases one has the dis- 
persionless spectrum with zero energy - the flat band 
in the vortex core and on the surface of the system 
correspondingly. 

This work is supported in part by the Academy of 
Finland, Centers of excellence program 2006-2011 and 
the European Research Council (Grant No. 240362- 
Heattronics). It is our pleasure to thank N.B. Kopnin 
for helpful discussions. 
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